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Recent preliminary results Irom the PAMELA experiment indicate the presence of an excess of 
cosmic ray positrons above 10 GeV. In this letter, we consider possibility that this signal is the 
result of dark matter annihilations taking place in the halo of the Milky Way. Rather than focusing 
on a specific particle physics model, we take a phenomenological approach and consider a variety 
of masses and two-body annihilation modes, including W^W' , , bb, t'^t~ , y,'^ , and e^e~. 

We also consider a range of diffusion parameters consistent with current cosmic ray data. We find 
that a significant upturn in the positron fraction above 10 GeV is compatible with a wide range 
of dark matter annihilation modes, although very large annihilation cross sections and/or boost 
factors arising from inhomogeneities in the local dark matter distribution are required to produce 
the observed intensity of the signal. We comment on constraints from gamma rays, synchrotron 
emission, and cosmic ray antiproton measurements. 



PACS numbers: 95.35.+d; 98.70.Sa; 96.50.S; 95.55.Vj 
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Dark matter in the form of a thermal relic is an appeal- 
ing explanation for the approximately 85% of the matter 
density of the universe not composed of baryons. In addi- 
tion to being a natural extension of the big bang cosmol- 
ogy, the candidates which naturally give the appropriate 
relic abundance have annihilation cross sections of the 
order of the electroweak scale, a natural scale for new 
particles in theoretical frameworks which provide a solu- 
tion to the hierarchy problem. 

If they exist, such thermal relics are expected to be 
annihilating in the halo today, generating potentially ob- 
servable fluxes of high energy particles, including gamma 
rays, electrons, positrons, and antiprotons. To this 
end, a number of cosmic ray and gamma ray experi- 
ments [ll, H, H, 0, 0, 0] have considered the search for 
dark matter annihilation products to be an important as- 
pect of their science mission. Of particular interest is the 
satellite-based cosmic ray experiment, PAMELA 0, Q. 
With its large acceptance (21.5 cm^sr) and excellent par- 
ticle identification, PAMELA is anticipated to measure 
the spectra of cosmic ray protons, antiprotons, electrons 
and positrons up to energies of 700 GeV, 190 GeV, 2 TeV 
and 270 GeV, respectively. 

Of particular interest for dark matter searches are 
high energy cosmic ray positrons [1, [l^, El and an- 
tiprotons [111 . \vX [Tsj . The spectra of such particles are 
generally expected to be dominated by the products of 
high energy cosmic ray interactions with the interstellar 
medium. In contrast, the spectra of protons and elec- 
trons are dominated by particles produced in astrophys- 
ical accelerators, e.g., supernovae. As a consequence, the 
ratios e"^/(e"^ -|- e~) and p/p are, in the absence of pri- 
mary sources of cosmic ray antimatter (such as nearby 
pulsars [13, [H, expected to fall at high energies. 

A signal of an upturn in these ratios would constitute 
strong evidence for a new primary source, such as dark 
matter annihilations. 



Early results of the PAMELA experiment [17[ show a 
dramatic upturn in the positron fraction from 10 to 50 
GeV (consistent with earlier indications from HEAT [l| 
and AMS-01 0), while showing no excess in the antipro- 
ton data. These data immediately invite interpretation 
within the context of dark matter [l^, [l^, [l^, [2l| . Such 
an interpretation, however, is not trivial. In particular, 
the observed positron spectrum is somewhat harder than 
the spectral shape expected from dark matter annihila- 
tions to hadronic modes. Furthermore, the amplitude of 
the signal is very large and potentially difficult to recon- 
cile with the expectations of a thermal relic. 

To calculate the cosmic ray spectra resulting from dark 
matter annihilations, we use the publicly available code, 
G ALPROP [22j . For a given choice of the diffusion coef- 
ficient, boundary conditions, energy loss rate, and cosmic 
ray injection spectrum and distribution, this code solves 
numerically the steady-state diffusion-loss eg nation (36j: 
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where dn/ dE is the number density of particles per unit 
energy, K{E) is the diffusion constant, and h{E, x) is the 
energy loss rate. The source term, Q{E,x), reflects the 
mass, annihilation cross section, dominant annihilation 
modes, and distribution of dark matter in the Galaxy. 

We have adopted three sets of diffusion parameters 
which we have found to provide provide good fits to 
B/C and sub-Fe/Fe cosmic ray data above 5 GeV, and 
Be^^/Be^ data above 1 GeV (for a review, see Ref. [23l|). 
Thoughout this letter, we will refer to the following as 
Models A, B and C: 

A: K{E) = 5.3 • lO^^ cmVs {E/4GcY)°-*^ , L =4 kpc 

B: K{E) = 1.4 ■ 10^8 cmVs (E/AGeVf-'^^, L =1 kpc 

C: K{E) = 7.3-1028 cmVs (E /AGcY)°-^^ , L ^6 kpc. 
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Model A 


Model B 


Model C 


Mass Mode 


xVd/ BF 


xVdf BF 


x'/df BF 


100 e+e" 
100 

100 r+r" 

100 w+w- 

100 ZZ 
100 66 


0.152 3.8 
1.028 6.1 
2.893 12 
1.758 24 
1.921 34 
5.154 33 


1.459 23 
0.175 25 
2.019 45 
0.728 91 
1.139 100 
4.692 100 


0.555 2.4 
1.577 4.3 
3.224 9.0 
2.259 17 
2.413 24 
5.107 24 


300 e+e" 
300 ^i+^i" 
300 T+T- 
300 W+W" 
300 

300 bb 


0.182 32 
0.186 44 
1.131 57 
2.598 66 
3.126 74 
4.133 57 


1.132 430 
0.475 250 
0.387 240 
2.483 240 
2.993 250 
3.735 180 


0.439 18 
0.532 29 
1.586 39 
2.781 47 
3.256 53 
4.216 42 


1000 e+e" 

lUUU ^ 

1000 T+r~ 
1000 W+H/- 
1000 ZZ 
1000 bb 


0.106 310 

0.333 430 
2.243 210 
2.552 210 
2.877 160 


1.533 6300 

0.118 2400 
1.757 820 
2.055 770 
2.270 570 


0.210 170 
u.ooy Z(U 
0.693 280 
2.515 150 
2.809 150 
3.141 110 



TABLE I: The quality of the spectral fit (x^ per degree of free- 
dom) and the boost factors required for various dark matter 
masses (in GeV), annihilation modes, and diffusion parame- 
ters to produce the PAMELA positron excess. The column 
BF contains the boost factors required assuming a local dark 
matter density of p = 0.35. As stated in the text, the X^/df 
should be interpreted as a qualitative distinction between the 
scenarios, as the data are still preliminary and errors only 
statistical. 



where L is the distance above and below the galactic 
plane at which charged particles freely escape the Galac- 
tic Magnetic Field. For the electron/positron energy loss 
rate, we adopt the GALPROP defaults which include 
losses due to synchrotron and Inverse Compton scatter- 
ing, with a local magnetic field of 5/iG. For the source 
term, we adopt an Narvarro-Frenk- White (NFW) halo 
profile [1^ with a local density of 0.35 GeV/cm'^. 

In Fig. [U we plot the positron fraction as a function 
of energy for various dark matter masses, annihilation 
modes, and diffusion parameters. In each case, we have 
normalized the dark matter annihilation rate to provide 
the best possible fit to the PAMELA data above 10 GeV 
(below which the effects of charge dependent solar 
modulation can be significant). Although the PAMELA 
data are still preliminary, publicly available interpreta- 
tions of these preliminary data have already appeared in 
the literature [H, [H, HO] , and we follow this approach. 
The normalization of the annihilation rate is proportional 
to the annihilation cross section, the square of the lo- 
cal dark matter density, and a quantity known as the 
boost factor, which parameterizes the effects of clumps 
and other inhomogeneities in the local dark matter distri- 
bution relative to a smooth NFW profile. In each frame 
of Fig. [Tl we adopt a default value for the dark matter 
annihilation cross section {av = 3 x 10~^^ cm'^/s) and the 
local density (0.35 GeV/cm'^), and vary the boost factor 
to obtain the required normalization. 



We find that dark matter annihilations to leptons (in 
particular to e~^e~ and /i"*"^") are naturally able to 
provide a good fit to the spectral shape observed by 
PAMELA. Annihilations to gauge bosons or quarks, how- 
ever, tend to produce too soft a spectrum. This can be 
ameliorated if the diffusion boundary is small enough to 
limit the contribution from more distant annihilations 
(diffusion model B). In Table I, we give a qualitative mea- 
sure of the quality of the fit to the PAMELA spectrum 
(the per degree of freedom for the data points above 10 
GeV) for each case. As the data do not include system- 
atic errors, and as they are presently highly preliminary, 
these numbers should not be taken as any rigorous mea- 
sure of acceptable fits, but simply a means to distinguish 
the different scenarios. We also show in Table I the boost 
factor required to normalize each case to the PAMELA 
spectrum, assuming a local dark matter density of 0.35 
GeV/em3. 

The annihilation rate required to generate the flux 
of positrons observed by PAMELA is quite large and 
thus can be constrained by cosmic ray antiproton, syn- 
chrotron [25I. \2^. and gamma ray measurements, espe- 
cially those of EGRET [23| . In particular, if the large an- 
nihilation rate is the result of the dark matter possessing 
a large annihilation cross section (such as a non-thermally 
produced wino in an anomaly mediated supersymmetric 
scenario, for example psi. [28j). then constraints on dif- 
fuse gamma rays from EGRET would likely be exceeded 
by a factor of two or more if those annihilations pro- 
ceed to heavy quarks or gauge bosons. This conclusion 
can be evaded, however, if annihilations proceed largely 
to leptons. If the halo profile is cuspy, annihilations to 
gauge bosons are tightly constrained by synchrotron ra- 
diation limits in the center of the galaxy [2^, [2^ . Both 
of these constraints can be evaded if the local dark mat- 
ter annihilation rate is boosted by inhomogeneities in the 
surrounding few kiloparsecs without boosting the anni- 
hilation rate throughout the remainder of the Galaxy. 

Although somewhat unlikely, it is possible that the So- 
lar System happens to be near a large dark matter sub- 
halo, leading to a large positron flux without the over- 
production of gamma rays or antiprotons throughout the 
halo [2^. Recent results from the Via Lactea II simu- 
lation found that although the overall annihilation rate 
throughout the halo is boosted by only a small value 
(^1.4) on average, there is small (~1%) chance that 
the local annihilation rate is enhanced by more than a 
factor of 10 as a result of a large nearby subhalo [30| . 
Recent results from the Aquarius Project are more pes- 
simistic, however [3l| . Alternatively, dark matter anni- 
hilating in a density spike surrounding a nearby black 
hole (M ~ lO^-lO'^Mg) could strongly boost the local 
annihilation rate [l^l- It is also interesting to note that 
the spectrum of positrons from a nearby subhalo or in- 
termediate mass black hole would appear harder than if 
the positrons were produced from throughout a larger 
volume, perhaps enabling a better flt to the PAMELA 
data for annihilations to bb and other non-leptonic chan- 
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FIG. 1: The positron fraction as a function of energy for various dark matter masses, annihilation modes and diffusion 
parameters, compared to the background from secondary production alone (bottom line). In each frame, the annihilation rate 
was chosen to produce the best fit to the PAMELA data above 10 GeV. The required boost factor was calculated using our 



default values for the annihilation cross section 



3 X 10 cm-^/s) and the local dark matter density (0.35 GeV/cm"^ 



nels m. 

In summary, the PAMELA excess of high energy 
positrons, confirming earlier excesses from HEAT and 
AMS-01, raises the exciting possibility that we are seeing 
evidence of dark matter annihilations. In this letter, we 
have considered a range of dark matter annihilation chan- 
nels and masses and find many scenarios which provide a 
good fit to the data. In particular, dark matter annihila- 
tions to leptons (especially e"'"e~ and quite easily 
fit the observed spectrum. Annihilations to heavy quarks 



or gauge bosons, in contrast, provide a poorer fit to the 
data. This can be improved if most of the annihilations 
occur locally (such as is expected if the Solar System re- 
sides near a large subhalo or if the Galactic Magnetic 
Field confines charged particles only to a region within 
1-2 kpc of the Galactic Plane). In almost every case 
we have considered, very large annihilation rates are re- 
quired to produce the observed signal. In particular, 100 
GeV (1 TeV) dark matter particles require annihilation 
rates boosted by a factor of approximately ~2.5 to 100 
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(~100 to a few thousand) relative to the rate expected 
for a typical thermal cross section {av « 3 x 10~^^ cvo? /s) 
and a smooth halo with a local density of 0.35 GeV/cm"^. 
Such boost factors could arise from inhomogeneities in 
the local dark matter distribution, such as the presence 
of a large nearby subhalo or a density spike surrounding 
a nearby intermediate mass black hole. 

Such a large annihilation rate could also be generated 
by a non-thermally produced dark matter candidate with 
a large annihilation cross section, but this scenario is 
tightly constrained by EGRET's measurements of the 
diffuse gamma ray spectrum (and to a somewhat lesser 
extent by cosmic ray antiproton and synchrotron mea- 
surements [IE nil)- For example, although non-thermal 
^100 GeV wino-like neutralinos could produce the mea- 
sured flux of positrons, they would also exceed the dif- 
fuse gamma ray constraint by at least a factor of two (see 
the discussion in Ref. [l^].) A non-thermal dark matter 
candidate which annihilates largely to leptons, however, 
could evade such constraints. 

Although one could argue that the current data dis- 
favor dark matter candidates which annihilate largely 



to quarks or gauge bosons, including neutralinos (for 
a possible exception, see Ref. [3), and prefer those 
which annihilate largely to leptons, including Kaluza- 
Klein dark matter in models with a universal extra di- 
mension [H,!!!] or "exciting" dark matter (XDM) [H], 
uncertainties in the diffusion model make such conclu- 
sions premature. As PAMELA data become available at 
higher energies, however, it will become increasingly pos- 
sible to discriminate between various dark matter models. 
Data from GLAST [23| will also be very useful in con- 
straining the possibility of annihilations to non-leptonic 
channels. 
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